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Spin dynamics of two-dimensional electron gas confined in an asymmetrical quantum well is stud- 
ied theoretically in the regime where the scattering frequency is comparable with the spin precession 
frequency due to the conduction band spin splitting. The spin polarization is shown to demonstrate 
quantum beats. If the spin splitting is determined by both bulk and structural asymmetry mecha- 
nisms the beats are damped at zero temperature even in the absence of a scattering. We calculate 
the decay of spin beats due to the thermal broadening of the electron distribution function and 
electron scattering. The magnetic field applied along the structure growth axis is shown to increase 
the frequency of the beats and shift system towards the collision dominated regime. 



I. INTRODUCTION 



The spin dynamics in semiconductor nanostructures 
is a field of topical interest [l[. The progress in the 
nanostructure technology allows to fabricate semiconduc- 
tor quantum wells (QWs) of high finesse where the carrier 
scattering time r caused by the disorder amounts to sev- 
eral tens of picoseconds. The conduction band spin split- 
ting A in such structures can reach several meV which 
implies At/S > 1. In this case the spin beats are ob- 
served in time-resolved experiments for low enough tem- 
peratures where electron scattering role is diminished Q ■ 



Theoretically, these spin beats were studied in Ref. [3j , 
where the consideration was restricted to the case of 
isotropic conduction spin splitting, see also Ref. [4k some 
results for anisotropic splitting were given in Ref. [5[ . The 
regime of isotropic spin splitting being comparable with 
the electron Fermi energy was discussed in Ref. [6] . How- 
ever, the evidence is growing that the spin splitting in 
the state-of-the-art QW samples is anisotropic function 
of the electron wavevector 0, Q . Therefore the full inter- 
pretation of the experimental data on low temperature 
spin dynamics in high-quality QW samples needs a the- 
ory which takes into account the anisotropy of the spin 
splitting. The goal of the present paper is to put forward 
the theory of spin beats for the case of arbitrary relation 
between the spin splitting and the collisional broadening 
and for the arbitrary anisotropy of the spin splitting. In 
Sec. [IT] we present the kinetic equation for spin distribu- 
tion function and discuss the dephasing of the beats due 
to thermal broadening of the electron distribution and 
due to the anisotropy of the spin splitting. Section [Hi] is 
devoted to the decay of spin beats due to the scattering, 
and the effect of an external magnetic field on the spin 
beats is discussed in Sec. HVl 
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II. KINETIC EQUATION AND SPIN 
DYNAMICS IN THE ABSENCE OF THE 
SCATTERING 

We consider a QW grown along z || [001] from zinc- 
blende lattice semiconductor. The spin dynamics of the 
two-dimensional electron gas in the absence of magnetic 
field is determined from the kinetic equation [j| [l(J, [ll], 
d 

-r— + s k xn k + Q k {s k } = 0, (1) 

at 

where s k is the average spin of the electron in the state 
with the wavevector fc, i.e. s k is the spin distribution 
function, tt k is the spin precession frequency originated 
from the fc-linear terms in the electron effective Hamil- 
tonian [M [H, Ell 

^so = t| (f -fife) = Pi{a y k y - a x k x ) + (3 2 (o-xk y - a y k x ) , 

with /?i and @2 being Bulk Inversion Asymmetry (BIA 
or Dresselhaus) and Stucture Inversion Asymmetry (SIA 
or Rashba) spin splitting constants , and Q k {s k } is the 
collision integral. Here we use the coordinate frame with 
x || [100], y || [010]. From now on we consider spin- 
independent isotropic scattering. We assume that the 
average spin splitting on the Fermi surface A = Kl kF 
(where hp is the Fermi wavevector, Cl kF is the angular 
average of \£l k \) is much smaller than electron Fermi en- 
ergy Ep . It allows us to avoid the antisymmetrization of 
the collision integral @, [l5[ and use the simple expression 

Qk{s k } = - [s k - s k ] , (3) 

T 

where s k = {2ir)~ 1 s k d(p k is the angular-averaged 
spin distribution function, (p k is the axial angle of the 



We ignore the k 3 terms in the effective Hamiltonian assuming 
that the carrier density is not too high. The Interface Inversion 
Asymmetry term (IIA) has the same symmetry as BIA term 
thus /3i can be considered as a sum of BIA and IIA spin splitting 
constants. 
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wavevector k and r _1 is the scattering rate. We note that 
the kinetic equation ([1]) itself is applicable for EpT 3> h. 

Our goal is to solve Eq. ([1]) for the arbitrary relation 
between spin precession frequency £lk and the scattering 
rate r _1 . We assume that at t = the quasi-equilibrium 
distribution of spin z-component is created 0, [l2| : 



s z ,k{t 



0) = -^{fh,x/2 ~ fk-1/2), 



where /fe )S = {exp [(fj, s — Ek)/ksT] + is the Fermi- 
Dirac distribution function, fi±i/2 are the chemical po- 
tentials for electrons with spin z-components ±1/2 re- 
spectively, Ek is electron's dispersion and k B T is the 
temperature measured in energy units. Under the as- 
sumption of low spin polarization, 

p _ Z)fe(/fc,l/2 ~ fk-1/2) 1 
Z/fe(/fe,l/2 + fk-1/2) 

the initial spin distribution reduces to 



«,,fc(0) 



k R T 



fk(l~fk), 



(4) 



where is the Fermi-Dirac distribution function with 
fi = (/X1/2 + fi-i/2)/2. For simplicity the temperature is 
taken to be zero, T = 0, the effects of non-zero tempera- 
ture will be briefly discussed below. 

First, we consider a limit of a clean system, where r — > 
00. Thus, we neglect the collision term in Eq. ([1]). Then, 
the spin of the electron with the wavevector k precesses 
around f2fc according to 



s*,fc(i) 



fe (0) cosfifci. 



(5) 



The total spin polarization can be represented as average 
of Eq. ([5]) over the initial spin distribution, which in the 
limit of T = can be replaced by an average over the 
axial angle ifk of the electron wavevector: 



: (0) J ^ cos (y\M + Jl - W1P2 sin 2 Vk t 

In the limiting case of the isotropic spin splitting where 
Pi = or fa = 0, i.e. either BIA or SIA term is present 
in the spectrum only, the absolute value of the vector n k 
is ifk independent and according to Eq. 

s z (t) = s z (0) cosil kF t, 

where kp is the Fermi wave vector. One can see that 
the total spin oscillates and returns exactly to the initial 
value each period 2-k/D,u f 0- 

The spin splitting anisotropy leads to a difference of 
the spin precession frequencies in different points of the 
k space. As the precession frequencies are not commen- 
surable the spin will never return to the initial value, 
however, its time decay will be very slow. For example, 
in the limiting case j3i = [3% = (which corresponds to 



the strongest spin splitting anisotropy) the angular in- 
tegration in Eq. ^ can be easily carried out with the 
result [f| 



!(t , = 8 ,(0,J„(^), 



(7) 



where Jo(a;) is the Bessel function. At (ikpt/h 3> 1 the 
spin oscillates and decays as i^ 1 / 2 . Such an asymptotics 
(i -1 / 2 ) holds true even in the case of the arbitrary ratio 
of Pi and P2 provided ^/ifhjhlt / fi ~> 1, the latter can be 
seen directly from Eq. ([6]) since the main contribution 
to the integral comes from the points of the stationary 
phase, ifk — 7r/2, 37r/2. We note that this relaxation is 
'reversible' as it is caused by the spread of the precession 
frequencies, like for the spin dephasing induced by the 
spread of g-factor values. 

In the end of this Section we shortly comment on the 
effect of non-zero temperature on the smearing out of 
the spin beats. For the simplest case of isotropic spin 
splitting the beats are damped at ksT Ep as 



«*(*) 



,(0) 



irat 



sinh (nat) 



(8) 



where a = flk F kBT/2Ep. On a long time-scale (at 1) 
the polarization oscillations decay exponentially with the 
time constant 



2E F 



(9) 



We note that the exponential time decay of the spin beats 
takes place in the regime of anisotropic spin splitting as 
well, the time constant has the same order of the magni- 
tude as Tp. If the initial spin distribution has an energy 
width AE > fcfiT the beats will be damped during the 
time ~ Ep /(rtk F AE), provided this time is shorter than 
the energy relaxation time. 



III. EFFECT OF SCATTERING 

Now we are in the position to discuss effect of the elec- 
tron scattering on the spin beats. First we solve Eq. ([1]) 
for the case of T = and afterwards we discuss the com- 
petition between the damping of the beats due to the 
scattering and due to the thermal broadening of the dis- 
tribution function. In order to solve Eq. ([T]) we rewrite 
it in components as 



ds x . k 

dt 
ds Vi k 

dt 



T 
T 



0, (10) 
0, (11) 



Qy,kSx,k ^x.kSy.k H - — 0- (^^) 



ds z . k 
dt 



In deriving Eqs. (|10p we made use of the fact that z is 
the principal axis of the system, i.e. the relaxation of the 
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spin parallel to z will not lead to the average in-plane 
spin polarization of the carriers. It follows from Eqs. 
(jTOf. (fTT|) that 



d 1 



d_ 

di 



d_ 



Sz,k 



— =0. 



(13) 

The initial conditions for Eq. ([13]) are s z (t = 0) = s z (0) 
and ds z (t = 0)/dt = which follows from Eq. (TH]) . 

In the isotropic case the spin precession frequency flk 
is angular independent and the last term of Eq. (fTe 
vanishes. The spin dynamics is then described by: 



According to Eq. (fT"4f one can identify two qualita- 
tively different regimes of the spin relaxation: (i) spin 
precession regime with the exponentially decaying spin 

beats at a frequency \J^\ F — (2/ T ) 2 and the decay time 

constant Tg — 2r and (ii) the collision dominated spin 
relaxation regime showing exponential damping of the to- 
tal spin with the spin relaxation time Tjyp — 1/ (&k F t) . 
The transition between the two regimes takes place at 
ttk F T = 1/2. 

In the general case, where both BIA and SIA terms are 
present, the spin precession frequency Qk is a function of 
(fk and all angular harmonics of s z ,fc are intermixed. Eq. 
(fTB")) can be reduced to a set of linear algebraic equations 
by the decompositon of s Zj k into angular harmonics as 
follows 



«»(i)=«,(0) 



cosh 



— sinh ( — 
q \2t 



(14) 



where q = y/l — 4f2 2 ,^T 2 , in agreement with Ref. 

Note that, in Ref. [6J] this expression [Eq. (17)] is pre- 
sented with an error: the second term in brackets is miss- 
ing and the condition ds z (t = 0)/dt = is violated. 



S Z ,fc 



E 



Eq. (|13|) is thus equivalent to an infinite system of the 
linear equations for the values s m (t) which reads 



J 



d_ 
di 



1\ ds z 



Of 



d_ 

di 



1 \ i - s n 



r 



i^2 s z,m+2 — 0. 



(15) 



Here the integer m runs from — oo to oo, S m fl is the Kro- 
necker symbol and 



H 2 



(Pi +0: 



Hi 



■0102 



(16) 



These equations should be supplemented with the initial 
condition s Zj7n (t = 0) = s z (0)<5 mj o and ds Ztrn (t = 0)/dt = 
0. 

Figure [1] presents the time dependence s z (t) obtained 
by the numerical solution of Eq. (fTS"]) . Different pan- 
els of Fig. Q] show the results for different ratios 02/01 
calculated at a fixed 0\. Eq. (fl6|) clearly shows that 
the dependence s z (t) is invariant under the replacement 
01 *-* 02- Different curves in each panel are calculated 
for different scattering rate t^ 1 . Even a small admixture 
of the other term to the strong first term in Hamilto- 
nian ©, e.g. 02/01 = 0.2, leads to the damping of the 
oscillations and quite complicated behavior of s z (t). 

The inclusion of the scattering first smears the oscil- 
lations. The damping time is non-monotonous function 
of r _1 in accordance with the qualitative discussion pre- 
sented above. For small scattering rates the scattering 
enhances spin relaxation, while for the stronger scatter- 
ing the spin dynamics becomes collision dominated and 
the spin relaxation time becomes longer. 



The possibility for the experimental observation of 
these oscillations depends strongly on the relation be- 
tween the average spin precession frequency tik F , scat- 
tering time r and thermal broadening time tt- Namely, 
the condition Qk F T > 1/2 guarantees the spin precession 
regime. Moreover, fifc F Tr should be larger than unity in 
order avoid thermal damping of the oscillations. 

The dominant broadening mechanism is therefore de- 
termined by the shortest of Tp and r. For very low 
temperatures tt is very long while r is determined by 
the carrier scattering on the interface roughness and re- 
mote impurities. In this regime the beats are damped 
due to electron elastic scattering. In the regime of the 
intermediate temperatures where electron-electron scat- 
tering dominates over the elastic processes governing the 
momentum relaxation [l6j , the scattering rate is propor- 
tional to t -1 oc k 2 B T 2 /HEp and scales as T 2 [12], whereas 
the thermal broadening rate r^ 1 oc ClkphpT / Ep oc T, 
see Eq. ©. Their ratio tt/t scales as kBT/h£lk F and 
can be both larger or smaller than 1. For example, in the 
experimental conditions of Ref. Q h&k F ~ 1 meV and for 
the temperatures larger than hQk F /kp ~ 12 K the damp- 
ing of the beats is due to the electron-electron scattering. 
We note that in the experiments quasi-equilibrium initial 
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FIG. 1: The time dynamics of the average ^-component of the electron spin for different ratio of BIA and SIA terms for a 
fixed /3i term. The time units are nh/ f3\kF- The ratio $2/ '/3i is shown above each panel. Various curves correspond to different 
scattering rates, see legend. 



spin distribution can be established during the energy 
relaxation time te which can be comparable with r. In 
such a case one can not separate two contributions to the 
beats decay. 

In the steady-state experiments one measures the spin 
polarization degree at the continuous spin generation 
rate. Let us assume that this rate G s is directed along pr- 
axis and is isotropic in the k space. Provided the carrier 
lifetime is much longer than the spin relaxation time, the 
stationary solution of Eqs. (fTUl) - ([T^|l reads 



; + G s T 



1 



(17) 



where the average spin-z component is given by 



Sz,k 



G s 



[(1 



n 2 r 2 



2n 2 2 r 2 )(i + n 2 r 2 - 2n 2 2 r 2 p - 1 



(18) 

We remind that f2p and f2| are defined by Eq. ([T5)| . For 
example, if the spin splitting is isotropic the proportion- 
ality coefficient between average spin s z and G s is sim- 
ply given by (f^r) -1 which coincides with Dyakonov- 
Perel' spin relaxation rate in the collision dominated 
regime (9. fl0|. 



IV. SPIN DYNAMICS IN A MAGNETIC FIELD 

It is well known that the magnetic field B || z slows 
down the spin relaxation in the collision dominated 



regime [17|, Il8j , The goal of this Section is to analyze 
the effects of the magnetic field for the arbitrary rela- 
tion between flk F and r . Here we neglect the non-zero 
temperature effects. 

The magnetic field has a two-fold effect on the spin 
dynamics. First, it causes the precession of electron spins 
due to Larmor effect and, second, the cyclotron motion 
of the electrons in the field leads to the rotation of the 
wavevector fc of a given electron and correspondingly the 
spin precession vector The effects can be taken into 
account as additional terms in the left hand side of the 
kinetic equation |T]), namely, x ojl and As& = u)c[k x 
ds /dk] with ui and loq being the Larmor and cyclotron 
frequencies respectively, see Ref. [18j. 

First, we start from the clean limit where the scattering 
is absent and the spin-splitting is isotropic. One can use 
the rotating frame of reference where the &k does not 
rotate. In this frame the electron spin experiences the 
magnetic field being a sum of three components: ttk, ijJl 
and effective field arising due to transition to the non- 
inertial frame which equals to ±up (upper sign for /?2 7^ 
and lower sign for /3j 7^ 0) The spin of the Fermi- 

surface electron exhibits the harmonic oscillations around 
the non-zero value oj 2 /(uj 2 + £l 2 ) with the amplitude 

/(co 2 



fl 2 ) and frequency 



where u> 



u>c iwi (sig n depends on the splitting mechanism in 
question) [201 ] - 

The difference of beats frequencies for BIA and SIA 
spin splittings in the magnetic field allows, in principle, 
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FIG. 2: The time dynamics of the average z-component of the electron spin in the magnetic field for different values of the 
scattering rate r _1 , see legend. Left panel: isotropic spin splitting P2/P1 = 0, right panel: strongest anisotropy P2/ Pi = 1. 



to determine experimentally the symmetry of the spin 
splitting. We note that, however, the cyclotron frequency 
is much larger than the Larmor one, loq 3> lol in the most 
conventional systems and we disregard the Larmor effect 
from now on. 

For the anisotropic spin splitting spin beats become 
anharmonic. Applying Eq. ^ to the case of strongest 
anisotropy (3\ = (3% = f3 and taking into account that k 
rotates with cyclotron frequency we obtain for the total 
spin 

«,(t) = .,(0)J (|^sbi^). (19) 

One can see that even a small magnetic field restores 
the strictly periodic spin dynamics with the period T — 
47r /u>c i-e. twice larger than the cyclotron one. The time 
dependence of the total spin is quite non-trivial: the spin 
beats demonstrate multiple harmonics. 

In the presence of the scattering, oscillations decay 
with time as it is seen from the left panel of Fig. [2] For 
the anisotropic spin splitting the situation is qualitatively 
the same, see the right panel of Fig. [2] 

If the magnetic field becomes so strong that loq 3> &k F 
the variations of the effective field fife occur on the 
shorter time-scale than the spin precession in this field. 
It means that the anisotropic part of the spin distribution 
Sfc - Sfc <C s Z) fe. In other words, spin rotation angles be- 
tween the strong variations of f2fc are small. This regime 
is analogous to the collision dominated regime where the 
scattering is faster than the spin rotation. Thus, Eq. 
([T]) can be solved by iterations and the total spin decays 
exponentially with the time-constant 

T s = (l+LU 2 c T 2 )/n 2 T, (20) 



where CIq is given by Eq. (fll)]) . We underline that Eq. 
(I2TJ1) (as well as the results of Ref. [18]) is valid for the ar- 
bitrary relation between Clk F and r _1 provided the mag- 
netic field is strong enough. 

The results presented above are derived for the case 
of the isotropic scattering, where all the angular har- 
monics of the distribution functions relax with the same 
time-constant r m — r (for m = 1,2, . . .). Qualitatively, 
these results hold true for the case of angular dependent 
scattering. Numerical solution of Eq. (fTJ) in the limit- 
ing case of the long-range classical potential, where the 
collision integral can be replaced by differential operator 
Qk{sk} = —T~ 1 d 2 Sk/dip1, shows almost no difference 
with the presented results. 

In conclusion, we have theoretically analyzed the spin 
dynamics of the two dimensional electron gas in asym- 
metrical QWs for the arbitrary relation between the scat- 
tering rate and the spin precession rate at low tempera- 
ture regime. We have demonstrated that the spin dynam- 
ics shows quite complicated beats for the anisotropic spin 
splitting. In this case the beats decay due to the spread 
of the spin precession frequencies. The suppression of the 
beats by the thermal smearing of the distribution func- 
tion, carrier scattering and the magnetic field has been 
studied. 
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